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We congratulate the authors on an excellent paper. It has inspired us to con¬ 
sider ways to incorporate QMC within SMC schemes in settings in which the 
transition density of the latent process is intractable and pseudo-marginal meth¬ 
ods are deploye d. In particular, c onside r filtering for pa r tially observed (jum p) 
diffusions (e.g. iFearnhead et all ( 20081) : Pollock ( 20131) : Pollock et all ( 20141) ). 
in which (in the simplest setting) the latent process is a diffusion satisfying the 
SDE, 


d Xt — ct(Xt) d t + dB t , Xq — x , t £ [0 , T]. 


( 1 ) 


In this setting the data comprises partial observations (arising at a finite col¬ 
lection of time points) of the latent process; the extension to noisy observation 
is trivial. The transition density_of t he latent p r ocess ( un der cer tain regularity 
conditions) can be shown ( Dachuna-Castelle and Florens-Zmiroul 19801) to have 
the following form between any two time points a and 6, where 0 < a < b < T, 
with Wpp b denoting the law of a Brownian bridge between x a and Xb over 
[a,b\: 


Pb-a{xb\x a ) =M(x b \ x a ,b-a) exp j J a(u) du 




OCpb-a(xb\x a ) 

( f b a 2 (W t )+a'(W t ) 
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( 2 ) 




To propagate particles between consecutive observation times, a and 6, one 
could simulate from the proposal p b -a(x b \xa.)-, perhaps by rejection sampling, 
and modify the weight of each particle by a factor corresponding to an unbiased 
estimate of (where W ~ W* a b ’ Xi> ). Supposing that 

vt S [0,T] (a 2 (Wt)+ a '(W t )) /2e[L,U], 


and letting k ~ Poi [{U — L) ■ (b — a)] and (£i,...) ~ U[a, 6], we have the rep¬ 
resentation: 


tp{W) = e - L ' (b ~ a) -E 
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An unbiased estimate of T(VF), using a finite dimensional realisation of the 

iid 


sample path, is obtained by sampling re and £1 


a simple Monte Carlo approximation (Beskos et al 


m _U[a, b] and employing 

2006). 




This scheme uses an unbiased estimator constructed by simulating re and 
then using a re-dimensional uniform random variable to approximate the inner 
expectation. Finding a lower variance unbiased estimator of ip(W) is desirable, 
and one would like to exploit RQMC. However, the dimension of the random 
variable being random, it is not straightforward to do this directly. One could 
instead sample re in the usual manner, and approximate the inner expectation 
conditionally using an RQMC point set. There is clearly a computational cost 
associated with such a RQMC method, which will only be appropriate for prob¬ 
lems in which the variance of the simple Monte Carlo estimator of is large 
and re is typically small. 
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